Noise-assisted classical adiabatic pumping in a symmetric periodic potential 
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We consider a classical overdamped Brownian particle moving in a symmetric periodic potential. 
We show that a net particle flow can be produced by adiabatically changing two external periodic 
potentials with a phase difference ip in time and x in space. The classical pumped current is found 
to be independent of the friction and to vanish both in the limit of low and high temperature. 
Below a critical temperature, adiabatic pumping appears to be more efficient than transport due to 
a constant external force. 



PACS numbers: Q5.40.-a, 05.40. Jc, 05.60.-k 
I. INTRODUCTION 

Recently there has been a considerable interest in small 
amplitude adiabatic pumping in mesoscopic electrical 
conductors Two oscillating, out-of- phase perturba- 
tions are applied which lead to small distortions of the 
shape of the system. As a consequence a directed current 
is generated Q . The pumped current is a consequence of 
quantum interference effects. An elegant formulation of 
the problem has been achieved by Brouwer ||] based on 
the modulation of the emissivities of the system Q . In- 
elastic scattering does not suppress the pumped current 
but introduces an additional more classical contribution 
to the pumped current due to rectification ||J^. For a 
broader view of this very active field we refer to a few 
recent works 

Adiabatic pumping is of general interest due to the fact 
that only very slow perturbations are required: further- 
more if the amplitudes are small the system remains at 
all times close to the stationary equilibrium state. Thus 
parametric adiabatic pumping can be viewed as a tool 
to investigate the near equilibrium properties of a sys- 
tem. Since perturbations can be applied locally such an 
investigation gives information on the system which can 
not be obtained through the application of global and 
stationary forces. 

It is the purpose of this work to complement the 
quantum mechanical discussions mentioned above and 
to investigate parametric pumping for a purely classi- 
cal system. We consider particles subject to damp- 
ing and thermal noise in a symmetric periodic potential 
Vo{x) = Vb(— s). In addition to the static potential Vo{x) 
two small amplitude time-dependent oscillatory poten- 
tials act on the particles. The perturbations we con- 
sider are periodic in time and are periodic in space with 
the same period as the static potential. We investigate 
the case where the perturbations have a double phase- 
difference both in time and space. As in the quantum 
case, a directed current is generated which is proportional 
to the frequency of the oscillating perturbations and pro- 
portional to the product of their amplitudes. A directed 
current results for almost all type of perturbations, un- 
less these perturbations have a special symmetry (phase- 



differences equal to a multiple of tt) . Interestingly, for the 
small amplitude perturbations considered here, the ther- 
mal noise is essential: the pumped current vanishes in 
the zero-temperature limit, is maximal at some interme- 
diate temperature, and vanishes in the high temperature 
limit. 

A symmetry breaking of the system is necessary to 
generate a directed current. In the noise assisted para- 
metric pumping process discussed here the symmetry is 
broken not with the help of the static potential Vq{x) but 
through the perturbations applied to the system. Non- 
equilibrium state-dependent noise with the same period 
as the potential but out of phase with a symmetric static 
potential also leads to a directed current [p2|-p6[. Thus 
this is an example of directed motion in a symmetric po- 
tential for which the symmetry is broken not by the static 
potential but only by the non-equilibrium noise. Simi- 
larly directed motion can be obtained in systems with a 
spatially symmetric potential but with a friction constant 
which is state dependent . Wc emphasize the sym- 

metry of the static problem, since typically, the recent 
literature has emphasized directed transport in systems 
in which already the static potential ps|-[2l|] is asymmet- 
ric Voix) 7^ Vo(— x). The examples discussed here and in 
Refs. p^-p^ demonstrate that static symmetric break- 
ing, i. c. the consideration of a rate/ief potential (22 -26 1 is 
not necessary, if either non-equilibrium noise, or pertur- 
bations applied to the system act in a symmetry breaking 
way. 

Quasi-adiabatic perturbations of particles subject to 
the Smoluchowski equation have recently been investi- 
gated by Parrondo . Below we present a discussion of 
small amplitude parametric pumping which closely fol- 
lows the discussion by Parrondo p7||. 



II. PARAMETRIC PUMPING 

The overdamped motion of a classical particle in an 
external potential and subjected to thermal noise is gov- 
erned by the Smoluchowski equation for the probability 
density p{x, t), 
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where p is the mobiUty, (3 the inverse temperature and 
J{x,t) — ~fiV' {x,T{t)) — jikTd / dx is the current opera- 
tor. Here the prime denotes derivative with respect to x. 
We consider a total potential V{x^T{t)) that is written 
as a sum of a symmetric periodic potential 

Vo{x) = Vb[l - cos(27rx/a)] (2) 

with period a, plus a perturbation 

AV{x,f{t)) = Xi{x)T^{t)+X2{x)T2{t). (3) 

Here Xi{x) and X2{x) are arbitrary spatial functions 
with period a and similarly Ti{t) and 12 (t) are arbi- 
trary functions of time with period l-Kjuj. We consider 
two special examples with purely harmonic driving. In 
both examples the time-dependent external perturbation 
is composed of two sinusoidal potentials with amplitude 
Va and Vb ■ In the first example the spatial functions act 
over the entire period, 

AT/(x,f(t)) = - VACOs(27rx/a)cos(wi) 

— Vb cos(27ra;/a -I- x) cos(ti;t -t- (ys) (4) 

with a phase difference x in space and a phase difference 
ip in time. In the second example driving is spatially 
localized at two arbitrary points x\ and xi in the interval 

[0,a], 

AV(a;, f{t)) ^ - Va cos{ojt)5{xi - 2TTx/a) 

~VBCOs{ijjt + (p)5{x2—'2,T:x/a) (5) 

where 5 is the Dirac delta function. In the following 
we assume that /^V{x,T{t)) changes slowly in time and 
that its amplitude is small compared to the unperturbed 
potential Va,Vb ^ Vq. 

The quantity of prime interest is the mean particle 
current, averaged over one period of space and one period 
of time. 



this limit, the system remains close to the adiabatic so- 
lution p'^{x,t), which is obtained by setting in Eq. (|l|) 
dp/dt = 0. The latter is given by an equilibrium Boltz- 
mann distribution 
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and therefore does not yield any current. (Here we have 
in addition to the adiabatic solution p^{x,t) introduced 
/oJ(a;,t) for later reference). Since the adiabatic solution 
is not associated with a current flow, we need to find the 
correction to this solution to determine the current. We 
expect that the correction to the adiabatic solution is of 
the order of the variation rate of the perturbation [ p7[ , 
which, for the potential Ay(a;,T(t)) we consider here, is 
given by u. We thus seek a solution of the Smoluchowski 
equation of the form 



p{x,t)^p^{x,t) + T{t) iJ{x,T{t)) 



(8) 



The correction T{t)v{x, T{t)), of order uj, to the adiabatic 
solution pQ^Xjt) gives rise to the non-vanishing particle 
current. Inserting the ansatz (^) into the Smoluchowski 
equation (|l|) and neglecting the time derivative of the 
correction, which is of the order w^, we arrive at 
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i7{x,f{t)) = V^(t)Po 



(9) 

This second-order partial differential equation for 
i7(a;, T{t)) has to be solved with periodic boundary condi- 
tions, £7(0, T(t)) = t7(a,T(t)), and the condition that the 
integral of v{x^T{t)) along the interval [0, a] vanishes. 
This second condition follows from the normalization of 
p{x,t) over one (spatial) period. We find 

X 
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Due to the periodicity of the potential in time and space, 
the second term of J(a;, t) does not contribute to the cur- 
rent. We begin by solving the Smoluchowski equation 
(m) in the limit of small driving frequencies, uj 0. In 
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where Ci (t) and C2 (t) are two vectors of integration con- 
stants. Explicitly, Ci{t) is given by 



Ciit) = J dxp+ix,t) J V^(,)Po (y.i) > (H) 




2 



where p^{x, t) is given by Eq. (^. The solution (^) of the 
Smoluchowski equation is then obtained by combining 
Eqs. (0) and (|0). Using the above solution, the mean 
current (||) can be easily calculated, 
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We now take the case AV{x,f{t)) = Ti{t)Xi{x) + 
T2{t)X2{x). With Vj,.,. = d/dTi+d/dT2 and Green's 



theorem, Eq. (12) can be rewritten as : 

A 
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dTi 
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This is the key result of our paper. Eq. (13) gives 
the pumped current in terms of the derivative of the 
adiabatic solution p^^Xji) and its companion pQ{x,t). 
It is valid for slowly time-varying periodic potentials 
AV{x,T{t)) of arbitrary shape and arbitrary strength. 
The expression Eq. (|l^) is most useful for the case of 
small amplitude perturbations. In this case the deriva- 
tives of p^{x,t) can be evaluated at zero amplitude and 
the integral over the area enclosed by the path / dTidT2 
is simply a multiplying factor. 



III. EXAMPLE WITH GLOBAL DRIVING 

We now specialize to the small amplitude regime. To 
do so, we can set Ti = r2 = in the integrand of 
Eq. (|l^). For the particular potential AV{x,f{t)) in- 
troduced in Eq. (^ the area enclosed by the pumping 
path is / dTi(ir2 VaVb sin sin X and a calculation leads 

A 

to the pumped current 
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2 PVo \l[^^y. LoVAVBsincpsinx , (14) 
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where Io{x) and Ii{x) are the hyperbolic Bessel func- 
tions of order zero and one, respectively. Equation (|lj) 
exhibits the main features of adiabatic pumping. We 
see that the adiabatically pumped current (/) is linear 
in the pumping frequency uj and the amplitudes Va and 
Vb of the two external potentials. The current is pro- 
portional to the sines of the temporal and spatial phase 
differences. An important consequence of Eq. ( p^ ) is that 
the current vanishes if either or % is a multiple of tt. 



This shows that a double phase difference, both in time 
and in space is necessary in order to rectify the noise. 
More generally, it can be shown that if the perturbation 
is written as product AV{x, f(t)) = X{x)T{t) with X{x) 
periodic in space and T{t) periodic in time, there is no 
pumped current, for any amplitude. An interesting sit- 
uation, which offers a simple physical interpretation, is 
the one for which the current is maximum. This happens 
when ip,x = ±f • By taking \Va\ = \Vb\, the perturba- 
tion can be rewritten in the form 



AV{x, T{t)) = Va cos(27rx/a ± ut), 



(15) 



the sign being determined by the relation between ip and 
X and between Va and Vb- The maximum current is 
hence generated by a traveling wave potential. It is to be 
expected that a traveling wave potential is a particularly 
efficient way of generating a current [^,^ . 

Let us now examine the temperature dependence of 
the particle current. It is given by the function 



/pump(" = PVo) = u 



h{u) 



(16) 



This function has been plotted in Fig.(P (solid line). We 
observe that /pump(u) vanishes both in the limit of low 
and high temperature and that it reaches a maximum at 
u ~ 1.426... This behavior can be understood as follows. 
At very low temperature (large u), thermal activation 
is negligible. Since the perturbation AV{x,T{t)) is fur- 
thermore very small, the particle remains trapped in the 
minima of the bare potential Vq{x) and there is no trans- 
port. For the small amplitude perturbations considered 
here this clearly demonstrates that there is no classical 
pumping without thermal noise. In the opposite limit of 
high temperature (small u), thermal fluctuations domi- 
nate and we have simple (symmetric) Brownian diffusion 
with zero average displacement. The maximum pumped 
current corresponds to a thermal energy of the order of 
the potential energy. 

It is further instructive to compare the adiabatic 
pumped current (|lj) with the current created by a small 
constant (time and space independent) external force F . 
In this case the overdamped Brownian particle experi- 
ences the potential, V{x) = Vq{x) — Fx. A similar cal- 
culation of the current, up to first order in F^ yields the 
following expression (see also [pO[), 



(I) 



fiF 



(17) 



The temperature dependence of this force induced cur- 
rent, given by /force(") = V^(")i is shown in Fig.(|l|) 
(dashed line). Here, in contrast to adiabatic pumping, 
the current is maximum at very high temperature. This 
is due to the presence of a small but finite slope in the 
potential. In the low temperature limit both currents are 
suppressed. 
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FIG. 1. Temperature dependence f{u — /3Vo), of the adi- 
abatically pumped current (solid line) and the current 
(^) generated by a constant external force F (dashed line). 

We notice, however, that for the force induced cur- 
rent, /force('") Me~^", the decay with temperature 
is asymptotically faster than for the pumped current, 
/pump(u) ~ e~^". As a consequence, there is a critical 
temperature /3c below which the pumped current is larger 
than the current generated by F. In this regime, adia- 
batic pumping is a more efficient transport mechanism 
than applying an external force. The value of the critical 
temperature Pc depends on both the area enclosed by the 
pumping path, the frequency and F. 



IV. EXAMPLE WITH LOCALIZED DRIVING 

In the example discussed above the perturbation po- 
tential is extended through the entire spatial period of 
the system. Is this a necessary condition, or is a spatially 
localized perturbation sufficient to generate a pumping 
current? To answer this question we consider a pertur- 
bation given by Eq. (H). As for the previous example, 
we can compute the current from Eq. (nSl) to obtain: 



(!) 



ujVaVb sin ip f /3Vb 



8^3 V^o^ 



h{f3Vo,x,,x2), (18) 



where h{f3Vo,xi,X2) = /(/3Vb, xi, 0:2) + /(-/3Vb, xi, 2:2), 
and 



f{PVo,Xi,X2) 



^(2/3Vo sin(a;i+a;2) sin(a;i — 0:2)) 



xi ~X2)+ / dxp„ (x) , (19) 



with Pq {x) — Pq{x,Ti = 0,T2 = 0). 6 is the Heavi- 
side step function. If we exchange xi and X2, we change 
the sign of the current, since it corresponds to a change 
of the sign of the temporal phase difference. The tem- 
perature dependence is similar to the one encountered 
previously (i.e. the current vanishes for both zero and 
infinite /3Vo). The interesting dependence of the current 



on xi and X2 is determined by /i(/3Vo = 1,2:1, 2:2) which 
is plotted in Fig.(^). No current is generated if xi = X2- 
The pumped current (the function /i(/3Vb, xi, ^2)) is dis- 
continuous along the line xi = X2- Of course this discon- 
tinuity disappears if instead of a delta functions slightly 
broadened functions are used to describe the perturba- 
tion. 

The reason that the localized perturbation leads to a 
pumped current is due to the normalization of the dis- 
tribution function. As a consequence even a localized 
perturbation generates a non-local response. 




FIG. 2. h{l,xi,X2), which describes the variation of the 
pumped current with respect to the positions xi and X2 of the 
delta perturbations. There is no pumped current for xi = X2. 

Fig. (^ gives an overview of the pumping currents that 
can be achieved with the two localized perturbations. To 
gain further insight it is useful to consider several cuts 
through Fig. (|). 

In Fig. (g|) (a) we keep the position of one perturbation 
at a fixed location X2 = (at the potential minimum) and 
consider the pumping current as a function of the position 
xi of the other perturbation. The pumped current is then 
maximal if xi is at the inflection point of the potential 
(xi — 7r/2). There is no current if the perturbation is 
located at the maximum of the potential. As xi increases 
further the current direction is reversed. As a;i passes 
through zero the current jumps from a negative value to a 
positive value. Note that in this case the pumped current 
is anti-symmetric around xi = t: since one perturbation 
is located at a symmetry point (2:2 = 0) of the potential. 
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(a) X2 = 



1 2 
Xl 




and is antisymmetric as a function of xi around this 
point. 

For global driving we found that the current is a si- 
nusoidal function of the spatial phase difference between 
the two perturbations. In contrast, in the local driving 
model considered here the pumped current depends not 
in a simple manner only on the spatial distance between 
the two perturbations but also on their absolute locations 
in the interval. This leads to the much more complicated 
behavior depicted Fig. (^) and Fig. (||) (b). 



V. DISCUSSION 




FIG. 3. Variation of the pumped current of the localized 
parametric pumping model. Shown is the function h for (a) 
X2 — 0, (b) X2 — ^ and (c) X2 — 2-k ~ xi. 

In Fig. (||) (b) we keep one perturbation at the inflec- 
tion point of the potential X2 = 7r/2. As a function of 
Xl the pumped current jumps at xi = 7r/2. The current 
decreases as xi increases past X2 = 7r/2 goes through a 
local minimum, reaches a local maximum for xi just be- 
fore Xl = TT. The maximal pumped current is achieved if 
Xl is a little to the right of the local potential minimum. 
Note that in this case we have no symmetry around 7r/2 
and neither are the two directions of current equivalent. 

In Fig. (^) (c) we consider the case where the two per- 
turbations are symmetrically located around the poten- 
tial maximum tt. The current then jumps at xi = X2 = 



We first point to the remarkable fact that the pumped 
current is independent of the mobility fi. This should be 
contrasted with the linear dependence on /i of a current 
generated with the help of a constant force. That the 
pumped current is indeed independent of fi can easily 
be seen by considering Eq. (^). The solution of this 
equation determines the correction i> of the density to 
the adiabatic solution: This correction is proportional to 
Since the current operator is proportional to fi the 
resulting net current is independent of /i. We emphasize 
that the pumped current is independent of fx not only for 
the particular models considered here but quite generally 
for all adiabatic perturbations. 

In summary, we have shown that a directed current 
can be generated in a symmetric periodic potential by 
adiabatically modulating two small external potential 
parameters. We have investigated a model with global 
(spatially periodic) perturbations and a model with lo- 
calized perturbations. For the model with global pertur- 
bations we find that a double — temporal and spatial 
— phase difference is necessary to generate a current. 
The maximum pumped current is obtained for a temper- 
ature kT of the order of the potential height Vq and for a 
perturbative potential corresponding to a traveling wave 
{(p,X = i'''/2). Similarly for the model with localized 
perturbations we find a pumped current unless the two 
perturbations are located at the same position or have 
a spatial difference of tt. Pumping arises through the 
subtle interplay between thermal fluctuations and cyclic 
variations of the potential. It therefore disappears in the 
limit of low and high temperature when either the po- 
tential or the thermal energy becomes predominant. We 
have demonstrated the existence of a potential dependent 
critical temperature below which adiabatic pumping is 
a more efficient than applying a small constant external 
force. The work presented here can be extended in difler- 
ent directions: Systems with open boundary conditions, 
several space dimensions and inertial effects are possible 
subjects for further research. 
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Nanoscale Dynamics, Coherence and Computation and 
by the Swiss National Science Foundation. 
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